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Abstract 

The electronic structure and magnetism of S^CelrOg, an Ir-based double perovskite system has 
been investigated using first-principles calculations. We found that a strong spin-orbit coupling 
dictate the electronic and magnetic properties of this system. A small value of U along with SOC 
could open up a gap at the Fermi level, offering the possibility of novel J e ff = ^ Mott quantum 
state. Our calculations reveal that the magnetic ground state is antiferromagnetic in agreement 
with the magnetization data and provide the value of spin and orbital moment for this system 
which is in agreement with the other isostructural Ir-based compound. 

PACS numbers: 71.20.-b, 71.70.Ej, 75.50.Ec 
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I. INTRODUCTION 



Until few years ago, the common belief has been that due to the extended nature of the 5d 
orbitals, electron correlation does not play any significant role in 5d transition metal oxides 
(TMO) and density functional theory (DFT) within local density approximation (LDA) or 
generalized gradient approximation (GGA) can explain the metallic ground state of these 
systems. However, recently there are reports of insulating ground state in some Ir-based 
oxides such as S^IrO^i, Na2irOjr, S^h^O^ etc. It has been shown that a combined effect 
of strong spin-orbit coupling (SOC) and onsite Coulomb correlation (U) can lead to an 
antiferromagnetic insulating ground state in such irridates with half-filled (5d 5 ) Ir 4+ ions^. 
Since the interplay of SOC and U dictate the physics of irridates, it has gained lot of attention 
in last few years and many exotic phases like quantum spin liquid, topological insulator have 
been realized in these systems^-. Double perovskite systems with formula A 2 BB'06 have 
attracted substantial attention during the last couple of decades due to many novel properties 
and it is therefore likely that a combined effect of SOC and U may lead to interesting 
properties in such systems. Recently within a model Hamiltonian approach, Chen et alM 
have shown the possibility of exotic phases induced by strong SOC in completely ordered 
double perovskite A2BB'06 where B' has either Ad 1 or 5d x electronic configuration. The 
importance of spin-orbit coupling in Ir based double perovskites La2CoIr06 and S^CoIrOe 
have been studied by Kolchinskaya et alM- using x-ray magnetic circular dichroism (XMCD) 
measurement and they find that the orbital moment is quite large and almost equal to the 
spin moment. In this context, Ir-based double perovskites may be interesting and warrant 
a detailed theoretical investigation. It would be more interesting to study a Ir-based double 
perovskite system where Ir is in 5d 5 electronic configuration which can then be a potential 
candidate for exhibiting novel J e //=| Mott ground state. The novel J e //=| Mott state has 
already been suggested in many other iridates with half-filled (5d 5 ) Ir 4+ ions like Sr 3 Ir 2 07^, 
CalrO^, S^IrO^ and this recently discovered quantum state can give rise to various kind of 
interesting low-energy Hamiltonians for instance the highly anisotropic Kitaev model which 
may be important for quantum computing^. With this motivation, we chose a well-ordered 
double perovskite S^CelrOg where Ir is in 4+ charge state (d 5 configuration) and Ir is the 
only magnetically active ion in this system. 

Sr 2 CeIr06 crystallize in monoclinic structure and shows an antiferromagnetic transition 
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FIG. 1: Crystal structure of S^CelrOe- 

at 21 K— . Neutron powder diffraction (NPD) study^ did not see any magnetic Bragg 
peak and hence concluded that the magnetic moment at the Ir site is very small. NPD 
measurement of isostructural compound La2CoIrOg^ where Ir is in d 5 state, also did not 
give any information of the magnetic ordering at the Ir site. However XMCD measurement^ 
on La2CoIr06 clearly see substantial magnetization at the Ir site and the estimated spin 
and orbital moments are respectively 0.20 /is and 0.18 \ib per Ir. This difference between 
NPD and XMCD results is attributed to the low neutron scattering cross section for Ir. 
Since neutron diffraction study did not provide any insight to the magnetic ordering and no 
XMCD measurements are available for Sr 2 CeIr06, a theoretical investigation will therefore 
be very useful to determine the type of magnetic ground state in Sr 2 CeIr06. In view of 
the above, we have investigated the magnetic properties of S^CelrOe within the density 
functional approach and find the system to be spin-orbit coupling driven Mott-Hubbard 
antiferro magnetic insulator where a novel J e // = | state can be realized. The remainder of 
the paper is organized as follows. In section [Til we shall describe the computational details 
and the crystal structure of S^CelrOe- Finally section UTTI will be devoted to the presentation 
and discussion of the results of our calculations. 
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FIG. 2: Partial Density of States for non-spin polarized configuration. 

II. COMPUTATIONAL DETAIL AND CRYSTAL STRUCTURE 

All the calculations reported in this work are carried out using full potential linearized 
augmented plane wave (FPLAPW) method as embodied in WIEN2k code^. The atomic 
sphere radii (.Rmt) of Sr, Ce, Ir and O are chosen to be 2.37 a.u., 2.22 a.u., 2.07 a.u., and 
1.69 a.u. respectively. To achieve energy convergence of eigenvalues, the wave functions in 
the interstitial region were expanded in plane waves with a cut-off RMTk ma x=7 , where Rmt 
denotes the smallest atomic sphere radius and k max represents the magnitude of the largest 
k vector in the plane wave expansion. The valence wave functions inside the spheres are 
expanded up to Z max =10 while the charge density was Fourier expanded up to G max =12. The 
Brillouin-zone integration was done with a modified tetrahedron method^ using 130 special 
k points in the irreducible part of the Brillouin-zone to achieve self-consistency. Exchange 
and correlation effects are treated within the DFT using LSDA with parameters of Moruzzi, 
Janak and Williams^ and a Hubbard U was included in the framework of LSDA+U^ 1 ^. 
A moderate value for the effective on-site Coulomb interaction XJ e ff — 2.0 eV is taken to 
describe the effect of correlation. The effect of spin-orbit coupling was treated using a 
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TABLE I: Energy of FM and AFM state along with the spin and orbital moments. Energy for FM 
state is assumed to be zero. For the LSDA+U+SOC calculations, orbitals moment of Ir is written 
in the parenthesis. 







Energy 
(meV/fu) 


Spin (Orb) Moment 
(Pb) 








Ir 


O 


Int 


Tot 


LSDA+U 


FM 
AFM 


0.0 
27.21 


0.57 
0.44 


0.06 
0.0 


0.21 
0.0 


2.0 
0.0 


LSDA+U+SOC 


FM 
AFM 


0.0 
-9.72 


0.44 (0.40) 
0.40 (0.42) 


0.04 
0.0 


0.21 
0.0 


1.53 
0.0 



second- variational scheme. 

Sr2CeIr06 crystallize with a monoclinic symmetry in the space group P2\/n. The de- 
viation from the ideal double perovskite structure is the tilting of the corner- sharing CeC>6 
and IrOe octahedra which break the cubic symmetry and also the elongation of the Ce-0 
bond length to 2.2 A at the expense of Ir-0 bond length which become 2.0 A. The tilting of 
the two octahedra in opposite direction make the Ir-O-Ce angle to be 156°. There is a little 
distortion in the octahedra in terms of angles but the bond-length with O ions are all equal. 
The crystal structure is shown in Fig. [TJ Unit cell contain two formula units. Coordinates 
of the two equivalent Ir, located at 2c site are (0.0, 0.5, 0.0) and (0.5, 0.0, 0.5). Since the 
distance between this two Ir is quite large (5.82 A), they could only interact via O-Ce-0 
path. All the calculations are performed with the experimental atomic positions and lattice 
parameters a = 5.8180 A, b = 5.8402 A, c = 8.2355 A, and = 90.225°^. 

III. RESULTS AND DISCUSSION 

We have analyzed the electronic structure of non-spin polarized LDA calculations to get 
an estimate of the charge state of Ce and Ir ions and the octahedral crystal field splitting. 
The calculated partial density of states (Fig. [2]) shows that both Ce-d and Ce-f derived 
states are completely empty. This confirm that in the ionic limit, electronic structure is 
consistent with the 4+ charge state of Ce in agreement with the experiment^ and hence Ce 
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FIG. 3: Partial Density of States of Ir-d of both the spin channels for (a) FM ground state within 
LSDA+U and (b) AFM ground state within LSDA+U+SOC. 

will be magnetically inactive in this system. Ir is therefore expected to be in 4+ (d 5 ) charge 
state. The partial density of states (PDOS) of Ir-d (Fig. [2]) reveal that the Ir-d e g (x 2 — y 2 
and 3z 2 — 1) states are empty while the ti g (xy, xz and yz) states are partially occupied 
consistent with the 4+ charge state of Ir and the crystal field splitting between t 2 g and e g 
state is appreciable (~2.5 eV). Further there is substantial hybridization between Ir-d and 
0-p states. Since IrOg octahedra are distorted and also tilted by an angle of 12° from the 
c-axis, this can lead to some admixture between t2 g and e g states. However Fig. |2] clearly 
reveal that the mixing between t2 9 and e g states are very little. Since t2 g and e g states are 
separated by a large gap (2.5 eV) and they hardly mix, one can assign an effective angular 
quantum number l e ff = 1 for the t2 9 manifold as in Sr 2 Ir04-. We also note that the tilting 
angle of IrOg octahedra in Sr 2 Ir04 is 11°—. 

We have next considered the magnetic properties of Sr 2 CeIr06. Our LDA calculations 
reveal that the magnetic solution can only be stabilized in the presence of Hubbard U at 
the Ir site. The results of our calculations for the ferromagnetic and antiferromagnetic 
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(a) LSDA+U 
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(b) LSDA+U+SOC 
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FIG. 4: Band dispersion along various high symmetry directions in the AFM state, calculated 
within (a) LSDA+U and (b) LSDA+U+SOC approach. 

arrangement of Ir spins are summarized in Table. [H From Table d we gather that the half 
metallic ferromagnetic state (see Fig. [3]- (a)) with integral magnetic moment is stable when 
the calculations are carried out in the framework of LSDA+U method. This can be easily 
understood as in the spin polarized calculation, the t2 9 states in the majority spin channel 
are completely occupied while it is partially filled in the minority spin channel leading to 
a stable half metallic state as illustrated in Fig|3]-(a). However upon inclusion of SOC not 
only the system become insulating (see Fig. [3]-(b)) but the AFM state become the stable 
solution in agreement with the experiment. Tabled reveal that the spin and orbital moment 
are almost equal in magnitude and parallel to each other, establishing the importance of 
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SOC as in other irridates. The insulating gap is calculated to be 30 meV. It is interesting 
to note that recent XMCD measurement^ for isostructural compound L^CoIrOg also see 
an equal spin and orbital moment along the same direction at the Ir site. 

In the following we argue that the importance of SOC and insulating solution can be 
rationalized by invoking J e //=| state. In Fig. 2] we show the LSDA+U and LSDA+U+SOC 
band dispersion for the AFM ground state. The LSDA+U band dispersion near Fermi level 
displayed in Fig. IH-(a) exhibit six bands from two Ir atoms in the unit cell, arising from the 
t2 g manifold. As expected for the d 5 configuration of Ir, in the majority spin channel the 
t 2g states of one Ir will be completely occupied while for the other it will be 2/3 occupied 
leading to a metallic AFM state. A significant change is observed in the band-dispersion 
upon inclusion of SOC (see Fig. IH-(b)) where we have a manifold of twelve SOC bands 
arising from the t2 9 up and dn states of two Ir. In the limit of strong SOC the t2 g orbitals 
with an effective quantum state l e ff = l form J e // = | quartet and J e // = | doublet. Here 
J e ff = \ doublet would be energetically higher than the J e // = | bands, since J e // = | is 
produced from the J5/2 manifold due to large crystal field splitting. Therefore we clearly see 
in Fig. 0]-(b) that a branch of total eight bands of two Ir which form J e ff = | manifold while 
the remaining four bands belong to J e ff = | manifold. Among the total 10 Ir-d electrons in 
the unit cell, 8 electrons completely occupy the J e // = § bands and hence J e // = \ bands are 
half-filled. Therefore a small but finite U can easily split these half-filled band into upper 
and lower Hubbard bands, yielding a Mott insulating state (see Fig. IU-(b)). We note that 
such a mechanism was suggested for S^IrO^ as well as CalrOy^ where J e // = | state was 
realized. We however note that the spin and orbital moments are calculated to be nearly 
equal in contrary to the expectation from a J e ff = \ state where the orbital moment is 
expected to be twice to the spin moment. The deviation may arise due to the tilting and 
the distortion of the octahedra as well as because of the Ir-d and O-p hybridization. 

In summary, our calculations find AFM ground state of S^CelrOe and establish that 
strong SOC dictate both the electronic and magnetic properties of this system. We also re- 
port the possible manifestation of a novel J e // = | Mott state in a Ir-based double perovskite 
system. 
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